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The mean time of contact formation between two ends of a protein chain shows power law dependence 
with respect to the number of residues, tcf ~ Fluorescence quenching measurements based on 
triplet-triplet energy transfer show variation in the value of scaling exponent a for different protein- 
solvent systems. This points to the relevance of the protein-solvent interactions (solvent quality) 
and hydrodynamic interactions in determining the time scale of contact formation. Here, starting 
from a non-Markovian diffusion equation supplemented with an exponential sink term that accounts 
for the energy transfer reaction between donor and acceptor groups, we calculate the mean time of 
contact formation using the Wilemski-Fixman closure approximation. The non-Markovian diffusion- 
reaction equation includes the effects of solvent quality and hydrodynamic interaction in a mean-field 
fashion. It shows that the contact formation dynamics is mainly governed by two time scales, the 
reciprocal of the intrinsic rate of quenching (fco"^)”^, and the relaxation time tq = pb^/ksT of the 
coarse-grained residue of an effective size b with solvent viscosity rj. In the limit of 1, the 

dominating effect of the reaction-controlled kinetics yields the scaling exponents as 0.89, 1.47 and 
1.79 in poor, theta and good solvents respectively. In the opposite limit k^'^ro 3> 1, the dominating 
influence of the diffusion-controlled kinetics results in a as 1.90, 2.17, 2.36 for a freely-draining 
and 1.31, 1.77, 2.06 for a non-freely-draining chain in poor, theta and good solvents respectively. 
In the intermediate limit, fco'^ro « 1, the increase in the number of residues switches the kinetics 
from reaction-controlled at low N to diffusion-controlled at large N. These general results suggest 
that experimental estimates of the scaling exponents reflect solvent-quality dependence of the mean 
contact formation time in the reaction-controlled limit. 


I. INTRODUCTION 
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The rate at which protein conformational space can be explored to form intramolecular contact between two residues in 
a polypeptide chain is an important elementary process in protein foldingii In the past few years, triplet-triplet energy 
transfer and photoinduced electron transfer between donor and acceptor groups located at two ends of a polypeptide 
chain have offered novel means to probe the role of intramolecular chain diffusion in determining the rate of contact 
formation and the time scale of protein conformational fluctuationsBy monitoring the triplet-triplet absorption, 
these experiments have obtained the effective rate constant kcF for contact formation between two ends of a protein 
chain, the reciprocal of which is equal to the mean time of contact formation, tcf = ~ 10 — 20, the latter 

follows a power law dependence with respect to the number of residues, tcf oc N°‘, where a is the scaling exponent. 
The typical values of a for different protein-solvent systems are 1.05 ± 0.06^ 1.5^*^^ and 1.7 ± O.li^ii^ for fV ~ 10 — 20. 
For relatively shorter polypeptides, where effects of stiffness become important, the dependence of tcf on N is much 

weaker 

In the asymptotic limit of large N, when the reaction between two end residues of a freely-draining chain in a theta 
solvent is considered instantaneous, the theoretical estimates of diffusion-controlled mean time of contact formation 
based on the Wilemski-Fixman (WF) and Szabo-Schulten-Schulten (SSS) formalisms yield oc and t^^^ oc 
respectively]The simplicity of the SSS formalism and its close agreement with the experimental scaling exponent of 
1.5 makes it a widely used theory to rationalize experimental data on end-to-end contact formation in polypeptides. The 
SSS prediction, however, deviates both from the WF theoryi^ and simulation resultsi^^r— In a recent work on contact 
formation kinetics^!, it is shown that the SSS formalism also yields t^^^^ cx N'^ once the monomer diffusion coefficient is 
replaced with an effective diffusion coefficient that includes the relaxation dynamics of the chain ends. The extended-SSS 
theory^! yields a ~ 2 in agreement with the WF formalism, but can not rationalize the weaker dependence {a < 2) of 
TCF on N observed experimentally. 

The WF formalism determines the mean time of contact formation by solving the reaction-diffusion equation in the 
presence of a sink term. The sink term accounts for the probability of end-to-end contact whenever the ends are within 
a contact distance a. For an idealized sink given by kj{R) = k^S{R — a), the reaction between two ends is instantaneous 
as soon as R = a. The diffusion-controlled (DC) limit of k^ —>■ oo, thus, yields ~ toN'^ as the mean time of 
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contact formation between two ends of a freely-draining ideal polymer/^ For measurements based on triplet-triplet 
energy transfer, however, the quenching rate depends exponentially on the distance between the donor and acceptor. 
The latter is given by kET{R) = exp(—2i?/a), where is the intrinsic quenching rate independent of the distance 
between donor and acceptor groups R and a is the contact distance for quenching. In the presence of a more realistic 
energy transfer sink, therefore, the reaction is not instantaneous but occurs at a rate that decays exponentially with R. 
Typical values of k^'^, range from 10® — For large but finite k^'^, it is conceivable that the weaker dependence 

of TCP on N can be rationalized using the reaction-controlled kinetics. Additionally, the presence of protein-solvent 
interactions which effectively amounts to change in the solvent quality from theta to good or poor solvent conditions, and 
hydrodynamic interaction which couples the dynamics of various residues in the chain can result in weaker dependence 
of TCP on N 

In this work, starting from a non-Markovian diffusion equation supplemented with an exponential sink term that accounts 
for the energy transfer reaction between two residues located at the exterior of a protein chain, we calculate the mean 
time of contact formation using the WF closure approximationThe non-Markovian diffusion equation describes the 
time evolution of the probability distribution of the distance between two residues on a protein chain and includes the 
effects of solvent quality and hydrodynamic interaction in a mean-field fashion. Our key result is that for triplet-triplet 
energy transfer, where the quenching rate depends exponentially on the distance between donor and acceptor groups, 
the kinetics of contact formation is reaction-controlled (RC) in the limit of <C 1 and diffusion-controlled (DC) in 

the opposite limit, kg^To :$> 1. Here, tq = rjb^/kpT is the relaxation time of the coarse-grained residue of an effective 
size b (molecular relaxation time). In the intermediate limit, k^^^To ~ 1, the increase in the number of residues switches 
the kinetics from reaction-controlled at small N to diffusion-controlled at large N. Our analysis shows that even for 
large values of the weaker dependence of top on N, can be rationalized using reaction-controlled kinetics in poor, 
theta and good solvent conditions. In the presence of the heaviside sink, the diffusion-controlled mean time of contact 
formation obtained from the present formalism is in agreement with the previous work based on generalized random 
walk description that accounts for non-local interactions approximately!^ 

This paper is organized as follows. Section II recapitulates the general features of a non-Markovian generalized Langevin 
equation (GLE), which is modified to include the effects of solvent quality and hydrodynamic interaction in a mean-field 
fashion. The GLE, when transformed into a diffusion equation and supplemented with an exponential sink term, results 
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in a non-Markovian diffusion-reaction equation. In Section III, the latter is used to determine the mean time of contact 
formation using the WF closure approximation. Section IV presents the main results of this calculation along with a 
brief discussion. Conclusions are presented in Section V. 

II. THEORETICAL MODEL 

The theoretical model presented below is based on a recent work where the time evolution of distance between donor and 
acceptor groups on a protein chain was described using an overdamped non-Markovian generalized Langevin equation 
(GLE) approach)^ This model was used to rationalize the results of a recent experiment where fluorescence quenching 
of photoinduced electron transfer between a pair of donor and acceptor groups on a protein chain was used to probe 
several universal aspects of protein conformational fluctuationsIn this experiment, distance fluctuations were shown 
to follow the Gaussian statistics with non-exponential decay, revealing non-Markovian nature of these fluctuations. The 
GLE approach^^ captured several universal aspects of the photoinduced electron transfer experiment, including correct 
prediction of the power law for the memory kernel^! and excellent agreement with two-point and four-point fluorescence 
correlation lifetimesi^ Here, we modify the GLE to include the effects of solvent quality and hydrodynamic interactions 
in a mean-field fashion: 


dt' Kmn{t - t')-^'R.rnn{t') = -KRm„(f) -f (1) 

dt' 

In the above equation Rmn = I’n — is the distance between two residues labelled as m and n and located at 
positions and r„ respectively on the protein chain. The first term on the right hand side is the effective elastic 
force due to chain connectivity with free energy F(Rmn) = where k = 3k bT/ (R^n) is the force constanti^ 

The variation in the solvent quality from theta to good and poor solvent conditions results in effective repulsive and 
attractive interactions between chain residues respectively. Within a mean-field description, the latter can be included 
by considering (R^„) = \n — where v is the Flory exponent with = 1/2, 3/5 and 1/3 for theta, good and 

poor solvent conditions respectivelyj ^^i^^ Thus the above GLE, while retaining the Gaussian and non-Markovian nature 
of distance fluctuations between two residues, accounts for the effective non-local interactions between chain residues in 



an approximate manner. 
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In Eq. ([1]), the mean and the variance of the Gaussian coloured noise fmn are given by (fmn(i)) = 0 and {fmn(t)fmn(^)) = 
ksT Kmn{t) respectively, where Kmn{t) is the friction kernel. The Laplace transform of Kmn{s) = /g°° dt exp{—st)Kmn{t) 
is given by 




^4^mn (s) 

1 ^4^mn (s) 


( 2 ) 


where (f>m.n{s) = /g°° dt exp{— st)(j)mn{t) is the dimensionless time correlation of distance fluctuations. The latter can be 
obtained from the Rouse and Zimm dynamics in theta, good and poor solvents. An outline of the derivation is presented 
in Appendix A. The final expression is 


(R„,„(t) ■ R,nn(0)) _ J2'^=i{cos{pTTn/N) - cos{pTTm/N)]‘^e 
(Rm„(0) • R™„(0)) J2'^=i{cos{pnn/N) - cos{pnm/N)]'^/p'^’'+^ 


where N is the total number of residues in the protein chain and Tp is the relaxation time of the pth mode.. The above 
equation accounts for the solvent quality and hydrodynamic interaction in a mean-field fashion. For freely-draining and 
non-freely-draining chain corresponding to the absence and presence of hydrodynamic interaction, the relaxation time 
is specified using and r^f, where superscripts R and Z refer to the Rouse and Zimm dynamics respectively. In terms 
of the longest relaxation time corresponding to the p = 1 mode, and Tp = Ti jp^'^, where r/* ~ 

and Ti ^ toN^" respectively. Here, tq rjb^/ksT is the relaxation time of the coarse-grained residue of size b (Kuhn 

length) with solvent viscosity 

Eq. O when transformed into a Smoluchowski equation is given by 


5P(R™„,<) 

dt 


ddrnnid^ 


dR„ 


R-nin PiR’nini^') T 


keT 


P{'Rmn,t) 


( 4 ) 


In an earlier work^ the survival probability of the unreacted donor state calculated from Eq. o yielded excellent 
agreement with another recent experiment^ measuring temporal decay of the transient absorption signals for fourteen 
mutants and wild type reaction center of protein dynamics modulated electron transfer reaction in early stage of 
photosynthesis. Starting from the GLE [Eq. (P)], a brief outline of the steps involved in deriving Eq. 0 are presented 
in Ref. (30). 
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Eq. (m is a non-Markovian diffusion equation, which describes the time evolution of the probability distribution of the 


distance between two residues on a protein chain with time dependent diffusion coefficient, Dmn{t) = — ■ In the 

presence of an energy transfer reaction between donor and acceptor groups at the exterior (or interior) of the protein 
chain, Eq. is supplemented with an energy transfer sink term, resulting in the following diffusion-reaction equation; 


5P(R„m,t) 

dt 




5R„ 


P(R 

mn 5 t) + 


keT 

« d^Ln 


P{R,nn,t) 


—kET^R mn )P(R mm t) 


(5) 


Triplet-triplet energy transfer follows Dexter electron exchange as a mechanism for fluorescence quenchingi^^ Thus, the 
distance dependent energy transfer rate expression is given by 


kET^Rmn^ — ^0 OXp ( ‘IRjyiji! u) , 


( 6 ) 


where k^'^ is the intrinsic rate constant which depends on the spectral overlap integral and a is the contact distance 
for quenching. In what follows, we use Eq. ([5]) to calculate the mean time of contact formation in the presence of the 
energy transfer reaction sink [Eq. (|6])] . 


III. THE MEAN TIME OF CONTACT FORMATION 

The mean time of contact formation can be obtained from the survival probability S{t) = f dRmnR(Rmn,t) that the 
donor and acceptor groups on the chain have not reacted at time t. From Eq. @, it follows that = — {k(t)), where 
{k{t)) = dRmn kET{Rmn) P{Rmmt). If S{t) is assumed to decay as a single exponential, S{t) Ri exp(—t/rcF) = 
exp(—kcFt), then the mean time of contact formation between two residues, which is the reciprocal of the effective rate 
constant for contact formation, is given by 


TCF 


k-^ - 



dtSit) 


(7) 









The solution of the diffusion-reaction equation, Eq. ([S]) can be written 
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P(R„„,t) = Peq{Rmn) “ / dt' 


dR-mn^i'^'n 




( 8 ) 


where G{Rmn,t\R'mn^^') conditional probability that the distance between two residues on a protein chain which 

was R'ran some initial time t' is Rmn at time t. Starting from Eq. dH), a closed form expression for G{Rmn, tlR^ru t') 
can be derived, the details of which are given in Ref. (30). Here, we simply state the final result: 


G (Rmn, i|RmnJ 0) — 


27r|n - 


3/2 


exp 


3 (R,nra — R(„ 


.it)r 


2\n - [1 - 


(9) 


From the above equation, the following equilibrium distribution can be obtained in the long time limit, t oo: 


7eg(Rmn) — 


27r|n — 


\ 3/2 

j exp 


3r: 


2 

'mn 


2\n — 


( 10 ) 


which correctly yields = \n — mp^6^ for the distance between two residues at the interior of the chain and 

for the end-to-end distance with n = N and m = 0. 

Eq. (|8]), when multiplied with kEriRmn) and integrated over dRmn yields 


pt pOO poo 

{k{t)) = (k)^^ - dt' dRrrm / dRm„ kET{Rmn) G{Rmn,t\R'^n,t') kET{R'mn) ^(R-mn:^0 

Jo J —oo J —oo 


( 11 ) 


where {k)^^ = dRmn kEriRmn) PeqiRmn)- Eq. (fTTI) is a nonlinear integral equation which can not be solved in 
a closed form. Using Wilemski-Fixman closure approximation we determine the mean time of contact formation, the 
derivation of which is discussed at length in Ref. (23). The final result is 


TCP 


k-^ 

f^CF 



((fc(0fc(0))e, - {k)lq) 

{k)lq 


= te + te, 


( 12 ) 











where 


{k{t)k{0))^^ = 


dRr, 


kET^Rmn) G'(Rmni ^ |R-mraJ 0) ^ET {Rmn) Req(Rrnn) y 


(13) 


is the time-correlation of the energy-transfer rates. In Eq. (nm, te and te are the reaction and diffusion-controlled 
contributions to the mean time of contact formation respectively. For end-to-end contact formation, with n = N and 
m = 0, the th and td are given by 


tr = {k}J 



/ 3 p 


dR exp(—2i?/a) g{R) 


(14) 


and 


Td 



(fc(0fc(0))eg 

{k)lq 


/r ir R e-2(«+^')/“ /(i?,t) 

3(/>(t) [1 - (/>2(t)]i/2 (J- dR e-2fl/“ g{R)f 


(15) 


where f{R,t) = exp 


3{R‘^ + R''^) 

sinh 

3RR'(p{t) 


3R^ 



•> 9\^) — 

2Af2^b2 


In the above equation, is 


given by 






(16) 


where and for the freely-draining and non-freely-draining chain respectively. In the next 

section, we numerically estimate tr and tr to obtain the mean time of contact formation using Eq. 



















9 


(a) Rouse Heaviside (b) Zimm Heaviside 




FIG. 1: Log-log plot for the variation of the diffusion-controlled mean time of contact formation td with N for a chain in good, 
theta and poor solvents. The results are obtained in the presence of the heaviside sink in the (a) absence (Rouse) and (b) presence 
(Zimm) of hydrodynamic interaction. Symbols are theoretical results and the solid lines are the corresponding linear fits showing 
the power law dependence of diffusion controlled mean time of contact formation with respect to the number of residues, td . 


IV. RESULTS AND DISCUSSION 


A. Diffusion-controlled mean time of contact formation with the heaviside sink 


In a previous theoretical study, the solvent quality dependence of the mean contact formation time has been obtained 
using a generalized random walk description which accounts for the non-local interactions in a freely-draining chain 
approximately!^ The nonlocal interactions have been included by modifying the connectivity term in the Edwards 
continuum representation of the polymer. This involves introducing a parameter h, with values 1/3, 1/2 and 3/5, which 
correspond to the average size of the chain in poor, theta and good solvents respectively. By solving the reaction-diffusion 
equation with the heaviside sink, the diffusion-controlled mean time of contact formation was shown to follow a power 
law scaling td oc . The third and fourth columns of Table 1 present the values of the scaling exponent (3 obtained in 
this previous study^l and in a recent Brownian dynamics simulations^! respectively. These values have been obtained 
in good, theta and poor solvents in the absence of the hydrodynamic interactions (freely-draining chain). 

To compare the scaling exponents obtained in the present work with the earlier work, we first obtain the diffusion- 
controlled mean time of contact formation in the presence of the heaviside sink in good, theta and poor solvents. This 
is done by replacing kET{R) in Eq. ([5]) with the heaviside sink given by knsiR) = ^hs for R < a and 0 otherwise. The 
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Solvent quality 

Present Work“ 

Previous Work^’“ 

Previous Work^’“ 

Present Work^ 

Good 

2.29 

2.28 

2.4 

2.0 

Theta 

2.08 

2.09 

2.0 

1.71 

Poor 

1.74 

1.76 

1.0 

1.2 


TABLE I: Comparison of solvent-quality dependence of scaling exponents f) for diffusion controlled mean time of contact formation 
td ~ in the presence of the Heaviside sink between present and previous^ work. Superscripts a and b refer to the solvent quality 
dependence of the freely-draining and non-freely-draining chain corresponding to the absence (Rouse) and presence (Zimm) of 
hydrodynamic interaction respectively. Superscripts 1 and 2 refer to previous estimates of the scaling exponents based on analytical 
theory and simulations for the freely-draining chain given in References [24] and [21] respectively. 


diffusion-controlled mean time of contact formation is calculated in the limit kns oo. In the absence and presence 
of hydrodynamic interaction corresponding to the freely-draining and non freely-draining cases respectively it shows a 
power law scaling oc which is depicted in Fig. (1). The values of the scaling exponents in good, theta and poor 
solvents in the absence of hydrodynamic interactions, obtained using this approach, are tabulated in the second column 
of Table 1. These values show excellent agreement with the previous analytical results [column 3]. The values of the 
scaling exponents obtained from simulations in good and poor solvents [column 4] are comparatively higher and lower 
than the previous and present analytical estimates. In Ref. (21), this has been attributed to finite size effects. 

The comparison of the diffusion-controlled mean time of contact formation obtained from the present and previous^! 
theoretical approach suggests that the non-Markovian diffusion equation approach, where the solvent quality is accounted 
for in a mean-field fashion, yields the same result as the generalized random walk description used earlier. However, the 
advantage of the present formalism is that it can easily be extended to account for the hydrodynamic interactions between 
chain residues. The scaling exponents for the diffusion-controlled mean time of contact formation in the presence of the 
hydrodynamic interactions are presented in column 5 of Table 1. In the presence of the hydrodynamic interaction, the 
scaling exponents are lower than the ones in the absence of this interaction [column 2]. This is because the presence of 
the hydrodynamic interaction couples the dynamics of different residues on the chain. As a result, the non freely-draining 
(Zimm) chain diffuses faster and has weaker dependence on N compared to the freely-draining (Rouse) chain. 







B. Mean time of contact formation with energy-transfer sink 
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The mean time of contact formation, in general, is a sum of reaction and diffusion-controlled parts given by Eqs. (HI 
and (fT51) respectively. While the integral in Eq. m can be evaluated analytically, there is no closed form analytical 
expression for Eq. (fTKll . We, thus, compute Eqs. (fTdl) and (ITKll numerically and substitute the result into Eq. (fT^ to 
obtain tcf- Before carrying out the time-integral in Eq. HSD, we non-dimensionlize time in terms of the time scale of 
intrinsic quenching by defining ti = tk^'^. 


Td= = 


poc 

/ dti 

{k(ti)k{0))^^ ^ 

/o 



(17) 


where is the dimensionless mean time of diffusion-controlled reaction. In terms of the dimensionless time 

ti, the modified expression for is given by 




(18) 


where k^ff = kg^^To is the dimensionless effective rate constant, tq is the relaxation time tq ^ rjb^/ksT of the coarse¬ 
grained residue of size b, j = 2v + 1 and 7 = 3^ for the Rouse and Zimm chain respectively. This non-dimensionalizes 
the time scale of contact formation = k^^TcF such as = Tp + Tp. The latter makes the integrals in Eqs. (IT4l) 
and p7p dependent on only through the dimensionless effective rate constant k^ff- 

Figs. (2)-(4) show the log-log plot for the variation of r^p with N in good, theta and poor solvents corresponding to 
V = 3/5, 1/2 and 1/3 for different values of fee//. For typical values of k^^ rs 10® — 10®s“^ considered in experiments, the 
dynamics of contact formation has a very sensitive dependence on solvent viscosity. For a fixed value of the solvent 
viscosity is varied by choosing different values of the dimensionless effective rate constant fee// = k^'^TQ. Figs. (2)-(4) 
show that the conditions kef / 1 and kef / ^ 1 correspond to the dominating influence of the reaction-controlled and 

diffusion-controlled kinetics respectively. In the former case when the solvent viscosity is low, k^'^ <C I/tq, the rate of 
molecular relaxation is faster than rate of reaction between donor and acceptor groups resulting in a reaction controlled 
kinetics. In the limit of high solvent viscosity, k^'^ I/tq, the contact formation kinetics is determined by the slower 


rate of molecular relaxation. 
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(a) Rouse Good 


(b) Zimm Good 




FIG. 2: Log-log plot for the variation of the dimensionless mean time of contact formation tqp with N for a chain in a good 
solvent {ly = 3/5). The diffusion-reaction equation [Eq. ([5ll] is solved in the presence of the exponential sink [Eq. ®] for the 
different values of the effective (dimensionless) rate constant k^ff = in the (a) absence (Rouse) and (b) presence (Zimm) 

of hydrodynamic interaction. Symbols are theoretical results and the solid lines are the corresponding linear hts showing the 
power law dependence of mean time of contact formation with respect to the number of residues, tqp ~ A^“. The limits of 
kf,f / 1 and fee// 1 show dominating influence of the reaction-controlled and diffusion-controlled kinetics respectively. In the 

intermediate limit, fee// ~ 1, the increase in the number of residues switches the kinetics from reaction-controlled at low N to 
diffusion-controlled at large N. 


In between these limits, Ki I/tq, the mean time of contact formation has significant contributions from both 
and r^. For small N, the end-to-end distance correlation relaxes to its equilibrium value fast, making it a reaction 
controlled process. At large N, the process is determined by the diffusion-controlled slow relaxation of the end-to-end 
distance correlation. The contact dynamics, thus, crosses over from the reaction-controlled kinetics at low N to the 
diffusion-controlled kinetics at large N. This is shown in Figs. (2)-(4) for good, theta and poor solvents in the absence 
(Rouse) and presence (Zimm) of hydrodynamic interaction. Lower the value of kgff, higher the value of N at which the 
crossover from the reaction-controlled to diffusion-controlled kinetics occurs. 

The values of the scaling exponents for the reaction and diffusion-controlled limits are indicated in Figs. (2)-(4) for good, 
theta and poor solvents respectively. The mean contact formation time for the reaction-controlled kinetics depends on 
the equilibrium probability distribution. As a result, the scaling exponents for the reaction-controlled kinetics in Figs. 
(2)-(4) are close to the values predicted by T^p ~ . The latter can be obtained by integrating Eq. (jlTl) analytically 

followed by an asymptotic limit of large N. The exact numerical integration in Eq. da, however, yield the values of 
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(a) Rouse Theta (b) Zimm Theta 




FIG. 3: Log-log plot for the variation of the dimensionless mean time of contact formation with N for a chain in a theta 
solvent {v = 1/2). The diffusion-reaction equation [Eq. ([5|] is solved in the presence of the exponential sink [Eq. ®1 for the 
different values of the effective (dimensionless) rate constant fce// = k^'^To in the (a) absence (Rouse) and (b) presence (Zimm) 
of hydrodynamic interaction. Symbols are theoretical results and the solid lines are the corresponding linear fits showing the 
power law dependence of mean time of contact formation with respect to the number of residues, tqp ~ A^“. The limits of 
kf,f / 1 and k^f / 1 show dominating influence of the reaction-controlled and diffusion-controlled kinetics respectively. In the 

intermediate limit, fee// ~ 1, the increase in the number of residues switches the kinetics from reaction-controlled at low N to 
diffusion-controlled at large N. 


scaling exponents which are slightly lower then Zv. 

In the limit of keff ^ 1, the dimensionless mean time of contact formation Tqp has dominating influence from the 
diffusion-controlled kinetics. In this limit. Figs. (2)-(4) show that the values of the scaling exponents in the presence of 
the exponential sink are slightly higher than the heaviside sink [Table 1]. This is because the scaling exponents for the 
heaviside sink have been obtained in the sole presence of the diffusion-controlled kinetics [Table 1]. For reaction-diffusion 
kinetics considered here, the limit of kgff ^ 1 ensures small but non-zero contribution from the reaction-controlled part 
with weaker dependence on N. The presence of latter modifies the values of slopes and intercepts in Figs. (2)-(4), 
resulting in slightly higher values of the scaling exponents. 

In the limit of keff « 1, as the solvent quality is varied from good [Fig. 2] or theta [Fig. 3] to poor [Fig. 4], the 
reaction-controlled kinetics occurs for a larger range of N. This is because in a poor solvent the difference in the mean 
contact formation time for the reaction-controlled and diffusion controlled limits is relatively smaller compared to theta 
and good solvent conditions. This requires higher values of N to attain the crossover. Similar trend is observed in the 
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(a) Rouse Poor (b) Zimin Poor 




FIG. 4: Log-log plot for the variation of the dimensionless mean time of contact formation tqp with N for a chain in a poor 
solvent {v = 1/3). The diffusion-reaction equation [Eq. ([5|] is solved in the presence of the exponential sink [Eq. ®1 for the 
different values of the effective (dimensionless) rate constant fce// = k^'^To in the (a) absence (Rouse) and (b) presence (Zimm) 
of hydrodynamic interaction. Symbols are theoretical results and the solid lines are the corresponding linear fits showing the 
power law dependence of mean time of contact formation with respect to the number of residues, tqp ~ A^“. The limits of 
kf,f / 1 and fee// 1 show dominating influence of the reaction-controlled and diffusion-controlled kinetics respectively. In the 

intermediate limit, fee// ~ 1, the increase in the number of residues switches the kinetics from reaction-controlled at low N to 
diffusion-controlled at large N. 


presence of the hydrodynamic interaction [Figs. (2b)-(4b)]. The reaction-controlled limit occurs for a higher range of N 
in the presence of hydrodynamic interaction than in its absence [Figs. (2a)-(4a)]. Again, this is because in the presence 
of the hydrodynamic interaction the diffusion-controlled mean time of contact formation are relatively smaller requiring 
higher values of N to attain the crossover from the reaction-controlled to diffusion-controlled limits. 


C. Comparison with experiments 

The end-to-end contact formation dynamics, as measured in different experiments, shows a power law dependence of the 
mean contact formation time with respect to the number of repeating units on a polypeptide chain with N ~ 10 — 20. For 
different polypeptide-solvent systems, the scaling relations have been found to be 1.05 ± 0.06,— 1.5^d^ and 1.7 ± 0.1— 
(good solvent). The scaling exponent of 1.5 has been rationalized using the SSS description of the diffusion-controlled 
contact formation dynamics in a theta solvent, which yields Tssb ~ £ The SSS prediction, however, deviates from 

several simulationsi^r— that predict Tqp ^ fV^, in consistent with the WF scaling—^ The mean time of contact formation 
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in the SSS theory is given by Tsss ~ /Dqu, where Dq = ksT/C^ is the monomer diffusion coefficient and C is the 

friction coefficient. In a recent study on diffusion-controlled contact formation dynamics, it has been shown that the 
SSS theory can yield the same result as the WF method if the monomer diffusion coefficient is replaced by the effective 
diffusion coefficient that accounts for end-to-end relaxation dynamics, thereby accounting for the higher-order modes of 
the chainji^i The latter yields the effective diffusion coefficient as He ~ resulting in the extended SSS prediction 

'T'esss in agreement with the WF prediction. This implies that the experimentally observed scaling exponent of 

1.5 can be rationalized using the SSS description provided the value of of Dq is considered to be much less than the 
monomer diffusion coefficientThus, the extended SSS prediction of Tesss ~ is in agreement with other theories 
and simulations, but can not rationalize tcf ^ scaling observed experimentally. 

While the above description in terms of the diffusion-controlled kinetics of contact formation assumes that the limit 
of ko ^ oo is always satisfied, the typical values of in different experiments lie between 10® — lO^s”^.'^ Based on 
the present analysis, the dominance of reaction or diffusion-controlled kinetics can be obtained by comparing 
with the time scale of molecular relaxation tq = rjh^/ksT. For typical values of rj ^ 1 — lOOcP^ T w SOOif and the 
coarse grained length of 5 2 — sA, an order of magnitude estimate yields fee// ~ 10“^ — 10“^ for k^'^ ss 10®s“^ and 

keff ~ 0.1 — 10 for 10®s“^. While the former (keff 1) corresponds to the limit where the reaction-controlled 

kinetics dominate, the latter {keff ~ 1) points to the crossover region where the reaction-controlled kinetics dominate 
at low N and the diffusion-controlled kinetics dominate at high N. The scaling exponents for the reaction-controlled 
kinetics obtained here are given by 0.89, 1.47 and 1.79 in poor, theta and good solvents respectively. The fact that they 
are close to the experimentally determined scaling exponents of 1.05 ± 0.06, 1.5 and 1.7 ± 0.1 points to the relevance 
of solvent quality dependence of the mean time of contact formation in the reaction-controlled limit. The importance 
of the reaction-controlled kinetics in rationalizing the values of the experimentally observed scaling exponents has been 
proposed earlieri^ 

Compared to the reaction-controlled kinetics, the dominating influence of the diffusion-controlled kinetics {keff ^ 1) 
yields a relatively stronger dependence of the mean contact formation on N. Figs. (2)-(4) show that in good, theta 
and poor solvents respectively, the freely-draining chains have mean time of contact formation governed by the scaling 
exponents 2.36, 2.17 and 1.90. For a non-freely-draining chain, in contrast, the respective values of the scaling exponents 
2.06, 1.77 and 1.31 in good, theta and poor solvents show comparatively weaker dependence on N. 


V. CONCLUSIONS 
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Starting from a non-Markovian reaction-diffusion equation that describes the time evolution of distance between two 
residues on a chain, we have calculated the mean time of end-to-end contact formation using the WF closure approxi¬ 
mation. This approach allows us to include the effects of solvent quality and hydrodynamic interaction in a mean-field 
fashion and yields a power law scaling of the mean contact formation with respect to the number of residues tcf ~ 

In the presence of the heaviside sink, the diffusion-controlled mean time of contact formation of a freely-draining chain 
in good, theta and poor solvents obtained from the present approach are in excellent agreement with the previous 
theoretical work based on a generalized random walk description. 

The non-Markovian reaction-diffusion equation when supplemented with a more realistic energy transfer sink shows that 
the interplay of reaction and diffusion-controlled kinetics determine the mean time of contact formation. In particular, 
the contact formation dynamics is governed by two time scales, the reciprocal of the intrinsic rate of quenching , 

and the relaxation time tq = r]b^/ksT of the coarse-grained residue of an effective size b. The limits of kg'^To 1 
and k^'^TQ ^ 1 show dominating influence of the reaction-controlled and diffusion-controlled kinetics respectively. In 
the intermediate limit, k^'^To w 1, the increase in the number of residues switches the kinetics from reaction-controlled 
at low N to diffusion-controlled at large N. These results show that the scaling exponent a has sensitive dependence 
on solvent quality mediated effective interaction between different residues on the protein chain, the solvent viscosity, 
the hydrodynamic interaction and the magnitude of the intrinsic quenching rate From these general results, we 

conclude that the experimental estimates of the scaling exponent a reflect solvent-quality dependence of the mean 


contact formation time in the reaction-controlled limit. 



Appendix A: Solvent quality dependence of the time correlation of distance fluctuations 
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The expression for the time-correlation of distance fluctuations in the absence and presence of hydrodynamics for a 
chain in a theta solvent has been derived in Ref. (28) using the Rouse and Zimm dynamics respectively. The latter also 
considers the case of a non-freely draining Zimm chain in a good solvent. Below, we extend this formalism to calculate 
the dimensionless time correlation of distance fluctuations [Eq. ([3])] for a freely draining (Rouse) and non-freely draining 
(Zimm) chain in good, theta and poor solvents. 

The dynamics of the nth residue at position r„(t) is described by the following Langevin equation^^: 


dvnit) 

dt 


Eh. 



(Al) 


where is the mobility matrix, U = Uei + Umt is the interaction potential which accounts for the entropic elasticity 
due to chain connectivity Uei = and solvent quality dependent effective interaction between 

chain residues, Umt- In a theta solvent Uint ~ 0. In addition, the absence of the hydrodynamic interaction implies that 
the mobility matrix takes the form The dynamics of a freely draining chain in a theta solvent is, thus, 

given by the Rouse dynamics 


dvnit) _ 3k bT d'^r„{t) ^ ^ 
dt ~ P dn^ ^ ’ 


(A2) 


where C, is the friction coefficient and fn(t) is the random force that follows the Gaussian statistics: (fn(0) = 0 and 
(f„(t) • im{t')) = 6kBTSnmS{t — t'). The condition of free chain ends require dr^/dn = 0 at n = 0 and n = N. Given 
this constraint, the position vector r„(t) in terms of the normal modes can be written as 

OO 

r^(t) = Xo(t)+ 2^Xp(i)cos(p7rn/A^), (A3) 

p^i 

where Aq = duVn. In terms of the normal modes, the Rouse equation is given by 


Cp 


dXpjt) 


dt 


kpXp{t) -I- fp(t). 


(A4) 








where Xp(t) = dn cos{pTTn/N)rn{t). The above equation can be easily solved to yield 
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(Xp(t) ■ X,(0)) = exp(-t/Tp). 


(A5) 


where Tp = Cp/kp is the relaxation time of the pth mode. In a theta solvent Tp = Ti/p^, where r/* = /Sir’^kET 

and kp = Gn’^ksTp^/Nb"^. The Rouse dynamics describes the dynamics of a freely-draining chain in a theta solvent. 
The effective attractive or repulsive interaction between chain residues in poor or good solvents respectively can be 
included in the parameter kp (in a mean field fashion) by determining kp = ■ The latter yields the solvent quality 

Wp/e« 

dependence of kp ~ p^~^'^''kBT/N‘^''b‘^ and 

In the presence of hydrodynamic interaction H„„ = ^ and Unm = —|- R„mR-nm + I 


for n ^ m. The nonlinear 


dependence of Hnm on Rnm makes the solution of Eq. (lAlIl difficult. The preaveraging approximation of Zimm replaces 
Hnm with {Hnm) Qq — J dR nm^^nm^eq ; wliGr6 Peq §iv6n by Eq. (j IQj) . This yi6lcis ~ iy (^|?2 tti^ Tigh'j — 

h{n — m)l. In terms of the normal modes Eq. (|Aip can be written as 


^ +f,(t)). (A6) 

9 

Following the linearization approximation^^, it can be shown that Cp = {hpp)~^ ~ rigbN'^p^~'^. Given that the solvent 
quality dependence of kp in the Rouse and Zimm dynamics is the same, the relaxation time of a non-freely draining Zimm 
chain is given by ~ jp^'^, where ~ N^'^b^rjg/kBT. This implies that the time correlation of (^p{t) ■ Xq(0)) in 

the presence of the hydrodynamic interaction obtained from Eq. (IMl) yield the same expression as Eq. with Tp 

replaced by . 

Since Rm„(f) = r„(t) — rm{t), it can be easily shown from Eq. (|A3p that 


4^mn (i) 


_(R 

mn (0-R rnn m 

(R rnn (0)-R mn m 


■ Xp(0)) [cos(p7rn/A^) - cos{pTrm/N)p 
(Xp(0) ■ Xp(0)) [cos(p7rn/7V) - cos(p7rTO/A^)]2 

^^i[cos(p7rn/A^) — cos(p7rm/IV)]^e“‘/'^r/pi+2i" 

^^i[cos(p7rn/iV) — cos{pTTm/N)]'^/p'^+'^'' 


(A7) 


which yields Eq. ©. 
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